Production of the exotic baryon in relativistic nuclear collisions 
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Recent experimental data appear to confirm the existence of the exotic <d + (uudds) baryon. Simple 
statistical considerations, which have generally proven rather successful, suggest that this particle 
should be produced in observable abundance in relativistic nuclear collisions at RHIC, where it 
may be identifed through the pK® invariant mass spectrum. The observation would not only add 
new evidence for its existence, but would also provide an additional means for probing the collision 
system, especially with regard to strangeness. 



PACS numbers: 25.75.-q, 12.39.Dc, 12.39.Mk, 14.20.-c 
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I. INTRODUCTION 

Based on an extended chiral soliton model, Diakonov 
et al. |]J made quantitative predictions for a previously 
conjectured exotic anti-decuplet of baryons, the lightest 
of which would be most easily identifiable. It can be 
considered as a uudds pentaquark resonance state hav- 
ing spin J = |, isospin 1 = 0, and strangeness S = +1 
(hence a charge of Q = +1). In the original paper 0, 
this particle was denoted Z + but Diakonov has more re- 
cently advocated that it be named instead. We there- 
fore employ this nomenclature. Its mass was prediced 
to be about 1530 MeV and its decay width less than 
15 MeV. The dominant decay modes are —> K + n and 
— > K°p. Various experiments have recently established 
support for the existence of this exotic baryon resonance, 
as briefly recalled below. 

Nakano et al. 0] used the Laser-Electron Photon Facil- 
ity at SPring-8 to irradiate 12 C with photons and consid- 
ered the reaction jn — > K~ K + n. They observed a sharp 
baryon resonance peak at 1540±10 MeV in the K~ miss- 
ing mass spectrum (corrected for Fermi motion). The 
Gaussian significance of the peak was reported as 4.6 a 
while the width was estimated to be less then 25 MeV. 
This evidence strongly suggests the existence of a mole- 
cular KN resonance with S = +1 that may be inter- 
preted as the 6. 

Barmin et al. employed a separated K + beam 
from the ITEP proton synchroton to irradiate the bub- 
ble chamber DIANA, filled with liquid zenon. They then 
considered the charge-exchange reaction K + n — > K°p on 
a bound neutron. This reaction can be fully measured 
due to the decay of the neutral kaon into a 7r~7r + pair 
inside the chamber. Their analysis of the K°p invariant- 
mass spectrum revealed a baryon resonance with a mass 
of 1539 ± 2 MeV and a width less than 9 MeV/c 2 . The 
statistical significance was estimated to be 4.4 a. This 
finding is a strong indication for the formation of the 0. 

Most recently, the CLASS Collaboration has reported 
results for photoproduction on deuterium and hydrogen 
at JLAB Analyses of the complementary missing 
mass spectra yielded a narrow S = +1 exotic resonance 
state in the nK + system (with a statistical significance 
of 5.2cr and 4.8c, respectively). The corresponding mass 



was found to be approximately 1540 MeV. 

In view of the mutually consistent experimental evi- 
dence for the existence of the 0, it may be reasonable 
to consider the implications for high-energy nuclear col- 
lisions where baryonic resonances are abundantly pro- 
duced [jj. Since the carries a unique combination of 
baryon number, charge, and strangeness it may provide 
a valuable additional diagnostic tool. 



II. STATISTICAL MODEL 

In order to help ascertain the practical prospects for 
observing the in relativistic nuclear collisions, we 
present here a rough estimate of its production yield. For 
this purpose we employ statistical considerations which 
have been quite successful in accounting for the observed 
yields of most hadronic species (see Refs. 0, H IE El ; 
for example). 

In a grand-canonical ensemble of hadron species i, the 
partition function factorizes into separate contributions 
for each specie, Z = TliZi, with 
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where plus is for bosons and minus is for fermions. The 
hadron mass is mi and <?i=2 Jj+1 is the spin degeneracy. 
The fugacity associated with the hadron specie i is 



Xi(T, { M }) = e ^u B+ Q^ Q+ S^ s )/T 



(2) 



where B i7 Q.j, and Si denote baryon number, electric 
charge, and strangeness, respectively. The ensemble is 
characterized by the temperature T and the three chem- 
ical potentials {/i} = {hb, A*q, Ms}- Finally, V is the 
enclosing volume of the system. 

Since the volume is hard to know a priori, it may be 
preferable to consider the corresponding spatial densities, 
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In the first line, the factors are arranged such that the 
terms in the sum are regular in the small-mass limit, 
mi — > 0, where K2 diverges. The second line exhibits 
the leading term, while the third line shows the simple 
exponential form that emerges when the mass is large, 
rrii ^> T. While we shall use the complete expression (01, 
it may be instructive to examine the quality of the three 
levels of approximation exhibited in Eqs. lolol) . 

Since the sign of the second term in J3J depends on 
the quantum-statistical nature of the particle, the lead- 
ing term overestimates fermion densities while it under- 
estimates boson densities. As it turns out, for all the pa- 
rameter values employed in our yield estimates, the first 
term in the expansion, Eq. Q, is an excellent approxi- 
mation for all baryon species, being accurate to within a 
few per mille. For kaons it is off by a few per cent, and 
even the pion densities are understimated by only around 
10%. The first term ||3J would thus be a quite reasonable 
approximation for rough estimates. 

By contrast, the large-mass approximation JSJ, which 
generally underestimates the densities, is never accurate 
for the situations of interest here. While underestimating 
the various baryon densities by less than «30% (only 
«20% for the 0), it is off by around half for kaons and 
underestimates the pions by about a factor of five. Thus 
it appears that this approximation is less useful. 



III. ESTIMATE OF THE YIELD 

We now employ the full expression for estimating 
the yield ratios that may be expected at RHIC. We fo- 
cus on the mid-rapidity region in central Au+Au colli- 
sions at the top beam energy of 100 GeV/A, where the 
highest freeze-out temperatures have been extracted. In 
this environment the net baryon, charge, and strangeness 
densities are quite small and, as one would expect, the 
extracted chemical potentials are relatively small. As we 
shall see, our estimates do not depend much on these 
parameters. In order to illustrate this insensitivity, we 
employ four different sets of chemical potentials. They 
are listed in Table H] and are motivated below. 

For the simpest estimate (A), we assume that all the 
chemical potentials vanish, fig = 0, [iq = 0, fj,s = 0. 
The second set (B) uses a finite baryon potential, fig — 
30 MeV, corresponding to values for which good fits have 
been obtained to the observed hadronic yield ratios [To| . 
(The sign of /is reflects the overall net baryon number.) 
In set C the charge potential has been given a finite value, 
fiQ = —4 MeV, which leads to g ood correspondance with 
the results quoted in Ref. [T(J]- (The negative sign of 
fiQ is a consequence of the initial neutron excess which 
introduces a corresponding isospin bias.) Finally, also 
guided by Ref. IQj, in set D the strangeness potential is 
finite as well, (13 = 10 MeV. (The fact that a non-zero 



value of /is is indicated by the data can be understood as 
a non-equilibrium feature reflecting the asymmetry in the 
elementary strangeness-producing reactions which causes 
kaons to be produced more readily than antikaons.) 



{/j,} set 


A 


B 


c 


D 


Ms (MeV) 





30 


30 


30 


fi Q (MeV) 








-4 


-4 


(is (MeV) 











10 



TABLE I: The four sets of chemical potentials employed in 
the calculations, = {/is, /j,q, /xs}. 

As for the temperature, we shall explore a broad range 
of values, T = 150 — 200 MeV, in order to bring out 
the dependence on this important parameter. (Most fits 
to the experimental yield ratios lead to chemical freeze- 
out temperatures around 175 MeV.) By comparing the 
results of the different parameter sets, we may then elu- 
cidate the sensitivity to the chemical potentials. 

We first calculate the ratio between the Q yield and 
those of certain familiar hadron species, namely A, p, 
K + , and tt + , which have different combinations of baryon 
number, charge, and strangeness. These yield ratios are 
listed in Table ITT1 for the four parameter sets A, B,C, D: 



T (MeV) 


150 160 170 180 190 200 


9/ A 

a e/ P 
e/K+ 

9/7T + 


9.0 10.7 12.5 14.3 16.1 18.0 
3.44 4.39 5.45 6.59 7.82 9.12 
0.74 1.12 1.62 2.23 2.98 3.85 
0.19 0.32 0.50 0.76 1.09 1.51 


6/ A 
B Q/p 
Q/K+ 

9/7T+ 


9.1 10.7 12.5 14.3 16.1 18.0 
3.44 4.39 5.45 6.60 7.82 9.12 
0.90 1.35 1.93 2.64 3.49 4.47 
0.22 0.38 0.60 0.90 1.28 1.75 


6/ A 
C Q/p 
Q/K+ 

9/7T+ 


8.5 10.2 11.9 13.7 15.5 17.3 
3.35 4.28 5.32 6.45 7.66 8.94 
0.90 1.35 1.93 2.64 3.49 4.47 
0.22 0.37 0.59 0.88 1.25 1.72 


9/ A 
D Q/p 
Q/K+ 

9/7T + 


9.1 10.8 12.6 14.4 16.3 18.1 
3.35 4.28 5.32 6.45 7.66 8.94 
0.90 1.35 1.93 2.64 3.48 4.47 
0.22 0.37 0.58 0.87 1.24 1.71 



TABLE II: Hadronic abundance ratios ne/rn (in per cent), 
with i = A,p, K + ,7r + , for a range of temperatures T (MeV), 
as obtained in the grand canoncial treatment ||2J with the four 
sets of chemical potentials {fi} listed in Tabic [I] 

It is evident that the calculated values show very little 
variation from one parameter set to the next, thus ren- 
dering the predicted yield ratios rather insensitive to the 
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precise values of the chemical potentials. Furthermore, 
taking the numbers at face value, the O yield should be 
about 12 — 14% of the A yield. Consequently, since good 
A data has been obtained, one would expect the 9 to be 
produced at a detectable level. 

In order to refine the estimate, one may further reduce 
the sensitivity to the chemical potentials by considering 
double abundance ratios for which the various attributes 
are counterbalanced. The leading fugacity factors then 
disappear (see Eq. (@J). This procedure may be particu- 
larly useful with regard to the strangeness which is not 
well described by the simple grand canonical model but 
requires that an additional suppression factor be intro- 
duced 0| ■ We therefore consider various double abun- 
dance ratios for which the values of the baryon number, 
charge, and strangeness balance out. The corresponding 
grand canonical results are listed in Table ITTT1 



M 


T (MeV) 


150 


160 


170 


180 


190 


200 


A 




13.7 


15.6 


17.5 


19.4 


21.3 


23.2 


B 


e/p 


13.7 


15.6 


17.5 


19.4 


21.3 


23.2 


C 


K+/7V + 


13.7 


15.6 


17.5 


19.4 


21.3 


23.3 


D 




13.8 


15.7 


17.6 


19.6 


21.5 


23.4 


A 




14.1 


16.0 


18.0 


19.9 


21.8 


23.7 


B 




14.1 


16.0 


18.0 


19.9 


21.8 


23.7 


C 


K°/w° 


14.1 


16.0 


18.0 


19.9 


21.8 


23.8 


D 




14.1 


16.0 


18.0 


19.9 


21.8 


23.8 


A 




8.8 


10.5 


12.3 


14.1 


15.9 


17.7 


B 


e/A 


8.8 


10.5 


12.3 


14.1 


15.9 


17.7 


C 


K+/K 


8.9 


10.5 


12.3 


14.1 


15.9 


17.8 


D 




8.8 


10.5 


12.3 


14.1 


15.9 


17.7 



TABLE III: Double hadronic abundance ratios (in per cent) 
for a range of temperatures T, as obtained in the grand canon- 
cial treatment with the four sets of chemical potentials {/j,} 
listed in Table Q 



For these double yield ratios, the dependence on the 
chemical potentials is practically negligible and the tem- 
perature dependence is considerably smaller than for the 
single ratios shown in Tabled thus making the "predic- 
tions" more robust. 

The first ratio shown, (n&/n p )/ {jik+ l n ir+)i is perhaps 
the most useful since it involves, in addition to the O, 
positively charged hadrons that can be directly detected. 
Since the experimental yields of those are well known, the 
tabulated numbers can be used to make absolute predic- 
tions for the O yield in terms of the observed yields of 
7r + , p, and K + . Inserting the multiplicities reported by 
BRAHMS 11], we then find, at mid rapidity for central 
Au+Au collisions at the top RHIC energy, 

dne ne/n p (dn K + / 'drj) dn p . 

dri n K +/n 7I + yan^+ldr^jQ dr\ 

i.e. about one per unit rapidity in each central event. 



The second ratio, {riQ /n p ) / (uk® /n^a), is included pri- 
marily for illustrative purposes. It has pratically the 
same value as the first ratio, as should be expected from 
the very small isospin dependence of the yields. But it 
is less readily observable (while the neutral kaons can be 
identified through their decay via K® to charged pion 
pairs, the neutral pions must be measured by means of 
their decay into photon pairs). 

The third reaction shown, (ne/n p )/(njf+ l n K°)i IS ac ~ 
tually easier to measure since p and K + can be detected 
directly, thus leaving only the task of identifying n^o 
through the — > 7r~7r + decay. (This identification is 
anyway required for the observation of the 0, sec below.) 
Alternatively, because of the very weak isospin depen- 
dence, one might fairly safely approximate n^-o by n^-. 

While the focus of the present brief note is on the top 
RHIC energy, where the yield is expected to be largest, 
the results given in Table III II may be useful for other 
scenarios as well since the double yield ratios depend 
primarily on the freeze-out temperature, with very lit- 
tle sensitivity to the various chemical potentials. 



IV. DISCUSSION 

As the above estimates suggest, the should be pro- 
duced with non-neglible rates in high-energy nuclear col- 
lisions. It may thus be worthwhile to extract its yield 
from the data. Obviously, the experimental observation 
of this exotic baryon in nuclear collisions would provide 
entirely independent evidence for its existence. More- 
over, if experimentally detectable, the would provide 
an additional means of probing the system, due to its 
unique combination of attributes. Especially useful may 
be the fact that it is a S = +1 baryon, which could help 
improve our understanding of strangeness production. 

One possible method of analysis for identifying the 
might be to consider the invariant-mass spectrum of the 
KgP system, as in Ref. 3|. In addition to the readily 
available kinematic information on the proton, this ap- 
proach requires the K® mesons to be identified through 
their ir~ir + decay vertex. This should be possible as 
well with the existing instrumentation (using the STAR 
TPC, for example). Presumably, since the intrinsic © 
decay width is so small, the observed width would pri- 
marily reflect the detector resolution (the construction 
of the p7r _ 7r + invariant-mass spectrum requires complete 
kinematical information for all three final hadrons). To 
the degree that the associated instrumental tolerances 
are well understood, it should then be possible to predict 
the observable width, thus facilitating the identification 
of the signal. 

The may also be of interest on the theoretical side 
since microscopic transport models usually attempt to 
include all known hadronic resonances. In these treat- 
ments, a primordial hadronic population generated by 
model-specific mechanisms undergoes further evolution 
involving both decays and reactions. As a result, both 
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the thermodynamic properties of the hadronic resonance 
gas (such as pressure and entropy) as well as the final 
yield of the detectable hadron species are affected by the 
selection of included resonances. 

Thus, insofar as the existence of the baryon reso- 
nance can be considered established, it would seem nat- 
ural that it be included in those treatments. While its 
expected relatively small abundance will likely render it 
ineffective with regard to modifying the overall thermo- 
dynamic conditions, it may play a role for the abundance 
of other relatively rare species. In particular, it might be 
of interest to investigate how strangeness-related quanti- 
ties might be affected by the presence of the 9. 

Just as this note was being finalized, it was hypoth- 
esized ^2 that the 9 may in fact not be an isoscalar 
but rather is an isotensor having five members. Three 
of them, 9 ++ , 9 + , and 9°, would have isospin violating 
strong decays, while the remaining members, 9 +++ and 
9~ would have weak decays and so be long-lived. This 
intriguing possibility gives the search for 9 production in 
nuclear collisions even larger urgency, since it could quite 
possibly settle this issue. Indeed, assuming that 9 ++ is 
produced at about the same rate as 9 + , as the statisti- 
cal model would suggest, it should be readily identifiable 
as a peak in the pK + invariant mass spectrum (a task 
simpler than identifying the 9 + ). 

V. SUMMARY 

Motivated by the recent apparent experimental confir- 
mation of the predicted exotic S = +1 baryon 9(1540), 
we have briefly discussed its relevance to high-energy nu- 
clear collisions. Obviously, if this particle indeed exists, 
it would enrich the hadronic dynamics and ought to be 
included in models treating the hadronic resonance stage 
of the collisions. 

In order to provide a rough estimate of the expected 
yield, we have invoked the framework of statistical equi- 
librium, which has been found to account very well for the 
observed hadronic abundancies over a wide range. These 
simple estimates suggest that the 9 may be produced at 



a level of approximately 6% relative to the proton yield 
and about 14% of the A yield, rather independently of 
the prevailing chemical potentials (as long as their values 
are reasonably consistent with the observed yield ratios). 

Further insensitivity to the somewhat poorly under- 
stood strangeness suppression effect can be obtained by 
taking double yield ratios designed to ensure counterbal- 
ancing of the various conserved attributes (baryon num- 
ber, charge, and strangeness). By this method, using the 
observed ir + , K + , and p multiplicities, we estimate that 
about one 9 per unit rapidity will be produced at mid 
rapidty in central Au+Au collisions at 100 GeV/A 

Since the expected 9 yield is thus not negligible, we 
encourage an analysis of the data for the purpose of es- 
tablishing 9 production in heavy- ion collisions. This task 
could be carried out by examining the invariant mass 
spectrum for the pR® system, with the K® being mea- 
sured by means of its ir~ir + decay vertex. In view of the 
possible existence of an entire 9 multiplet [lj, ^ might 
be worthwhile to examine also the pK + invariant-mass 
spectrum at the same time. 

The observation of such a resonance in high-energy nu- 
clear collisions would provide independent supportive ev- 
idence for the existence of this exotic particle. Once its 
existence has been established, the 9 may offer an inter- 
esting additional means for probing the properties of the 
collision system (especially those related to strangeness) 
due to its unique combination of baryon number, charge, 
and strangeness. 
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